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ABSTRACT 

In many cosmological models, the large angular scale anisotropy in the 
cosmic microwave background is parameterized by a spectral index, n, and a 
quadrupolar amplitude, Q. For a Peebles-Harrison-Zel'dovich spectrum, n — 1. 
Using data from the Far Infra-Red Survey (FIRS) and a new statistical measure, 
a contour plot of the likelihood for cosmological models for which — 1 < n < 3 
and < Q < 50 fiK is obtained. We find that the likelihood is maximum at 
(n,Q) = (1.0, 19 fiK). For constant n the likelihood falls to half its maximum 
at Q pa 14 /iK and 25 fiK and for constant Q the likelihood falls to half its 
maximum at n ~ 0.5 and 1.4. Regardless of Q, the likelihood is always less than 
half its maximum for n < —0.4 and for n > 2.2, as it is for Q < 8 /iK and 
Q > 44 /iK. 
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1. Introduction 



With the COBE/DMR detection of anisotropy in the cosmic microwave background 
radiation (CMB) flBennett et al. 1992| ; ^moot et aTTggg ; |Wright et al."T9^ ), followed 



by the FIRS-COBE/DMR cross-correlation flGanga et al. 1993|) and a host of recently 



reported detections on smaller angular scales, we must now characterize the nature of the 
detected variations. 

The inflationary scenario predicts a spectrum of primordial density perturbations of 
form P(k) oc k for large angular scales ( fSteinhardt fc 'turner 198~3] ), while popular variants 
of cold dark matter (CDM) models predict P{k) oc k n with 0.7 < n < 1 (see Turner (1993) 
for a review). Textures will also appear scale-invariant (i.e., n = 1; [Fen et al. 1993j ), while 



cosmic strings may mimic an n of 1.35-1.5 ([Perivolaropoulos 1993 ). Thus, by limiting the 



range of n allowed by the data we can constrain theories of structure formation. It should be 
pointed out, however, that many theories are not well parameterized this way {e.g., Ratra & 
Peebles (1994)) and that in this work we consider models with Gaussian fluctuations only. 

If one assumes this power law spectrum for density perturbations in an Einstein- 
DeSitter universe, then the angular power spectrum of anisotropies in the CMB is given by 
( Pond fc Efstathiou 1987)) 



Q = Q 



2 47rr((2/ + n-l)/2)r((9-7i)/2) 
5 r((2Z + 5-n)/2)r((3 + 7i)/2))' 



where Q = ^J5C2/4:TC is a parameter that sets the amplitude of the theoretical angular 
power spectrum and correlation function at the angular scale of the quadrupole. [| Note 
also that Equation [I] does not account for microphysical processes that can change the 
power spectrum of anisotropies and thus alter the effective n; for example, "standard" scale 
invariant CDM has an effective n of 1.15 for /<30 ( Pond 1994Q . 



The first estimate of n using large scale CMB anisotropy data was made by fitting 
the theoretical auto-correlation function (ACF) derived from Equation [I] directly to the 
auto-correlation function of the COBE/DMR first year data. The result was n = 1.15l ;65 
and Q = 16.3 ± 4.6 /iK ( [Smoot et al. 1992j ). Seljak and Bertschinger (1993), using 



5 Smoot et al. (1992) use the notation Qrms-ps to denote the parameter, but for the sake of brevity we 
follow the notation of Wright et al. (1994b) and use Q. Note that Q is not the root-mean-squared amplitude 
of the quadrupole in our Universe; we do not measure the intrinsic quadrupole because of our limited sky 
coverage. 
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a maximum-likelihood analysis again on the first-year COBE/DMR, data, found that 
Q = (15.7 ± 2.6)e 46 ( 1 - n ) /iK. Additionally, Smoot et al. (1994) and Torres (1994) have 
both done topological analyses of these data to limit n to 1.7lo!e an d 1-2 ± 0.3, respectively. 
Bond (1994) has extracted the index from an analysis of the power spectrum of both FIRS 
and DMR. He finds n = 1.8±g;| for FIRS and n = 2.0±g:| for the first-year DMR maps. 
Wright et al. (1994b), using a power spectrum analysis with the two-year DMR maps, 
find n = l-46i ;44- Bennett et al. (1994), also using the DMR two-year maps, perform a 
maximum-likelihood analysis with Monte Carlo simulations and find n = 1.59io; 55 . Finally, 
Gorski et al. (1994) use a power spectrum based method to compute the likelihood of 
models. Applying it to the DMR two-year data, they obtain a maximum likelihood at 
(n, Q) = (1.2, 17 fiK) and a maximum in the marginal likelihood for n at 1.1 ± 0.3. 

This range of results, often with the same data, indicates the difficulties involved in the 
analyses. There is clearly a need for more data and consistent statistical techniques. In this 
Letter, we describe a method for analyzing CMB anisotropy maps and apply it to the Far 
Infra- Red Survey (FIRS). This method uses the correlation function of anisotropy maps. 
Thus, it is relatively insensitive to noise contamination that can infect analyses based upon 
the x 2 or the power spectrum of a map ( Bond 1994j) . 



2. The Far Infra- Red Survey 



FIRS is a balloon-borne, bolometric anisotropy experiment that observes in four 
channels at 170, 290, 500 and 680 GHz. It has a beam full- width-at-half- maximum 
(FWHM) of 3.°8 and has been flown successfully twice, resulting in coverage over most of 
the Northern hemisphere. The data used here are from the 170 GHz channel of the October, 
1989 flight and cover roughly one quarter of the sky. The data are shown in Ganga et al. 
(1993). The experiment and map are described more fully in Meyer et al. (1991) and Page 
et al. (1989, 1990, 1992, 1994a,b). 

For the purposes of this Letter, it is sufficient to note that the map is a set of pixels i, 
each with an associated residual temperature ti and statistical weight W{. The t^s are the 
average temperatures of all measurements made within pixel i after the best fit offset, dipole 
and Galactic model (the IRAS 100 /jm map smeared to a 3.°8 beam width) are removed. 
The statistical weights are the reciprocals of the variances of the measurements for each 
pixel. A 15° Galactic latitude cut has been imposed on these data and on the simulations 
described below to minimize the effects of residual dust emission near the Galactic plane. 
For this analysis, we assume that all the correlations in the map are due to correlations 
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in the CMB. In Page et al. (1994b), the data reduction and known systematic errors are 
discussed. There are no known systematic effects that introduce correlations into the data 
at levels which affect the results of this analysis. 



3. Analysis 

The correlation function of the CMB is estimated from the data with 

C(9 A ) = -i- J2 WiWjUtj, (2) 
A ijeA 

where the sum is over all pixels i and j for which the angle between them is 6a and where 

T A = E w i w r ( 3 ) 

C(6a) is presented in Figure 1. The correlation function is divided into 64 2.°8 bins, of 
which, because these data do not cover the entire celestial sphere, only n c = 58 contain 
data. Also, the correlation function at 6 = (or, alternately, the rms of the data) is 
excluded from the analysis, as it is affected most by instrumental noise. 

Simulations of the sky are made for each n from 0.9 to 2.9 in steps of 0.1 following 
the methods of Cottingham (1987) and Boughn et al. (1993). Specifically, maps of unit 
amplitude are made with underlying temperatures at each pixel i defined by 

T *= E E a^W^M,^), (4) 

1=2 m=-l 

where W\ is related to the Legendre transform of the experiment's beam response. If one 
assumes that the response is Gaussian, as we do here (but see Wright et al. (1994a)), this 
window function can be approximated by W\ = e ' where 61/2 is the FWHM 

of the beam (in radians). With a FWHM of 3.°8, Wi falls to 0.5 at I = 29. For these 
simulations, l max is either 150 or the highest / for which W1C1/W2C2 > 10~ 7 , whichever is 
lowest. Each a^o is drawn from a normal distribution with a variance of Ci/(0.8ttQ 2 ) and for 
m > both the real and imaginary parts of a^ m are drawn from normal distributions with 
variances of Ci/(1.QttQ 2 ). Finally, ai- m = (—l) m a* l m . We call these signal maps. In order 
to include experimental uncertainties, separate noise maps are created. The temperature 
at each pixel i is drawn from a normal distribution with a variance of The pixels in 

both the signal and noise maps have the same weights as the real data. In order to match 
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the processing done on the real data, an offset and a dipole are removed from each map 
separately in a least-squares fit. 

Note that if the maps contain nothing but statistical noise, Ta is the reciprocal of the 
variance in bin A. That is, o~c(e A ) = V^a (see Ganga et al. (1993) for a derivation and 
Smoot et al. (1994) for a discussion of weighting schemes for different statistics over CMB 
maps). 

With this prescription, we form C S s(9), the ACF of a signal map, C S n(9), the 
cross-correlation between a signal map and its associated noise map, and Cnn(0), the ACF 
of a noise map. We note that the full ACF of a single sky realization with a particular value 
of Q and noise will be 

4-7T / ATT 

C{6) = Q 2 —C ss (9) + 2Q^—C SN {6) + C NN {6). (5) 

This scaling allows one to reduce the number of simulations made, though a set of 
simulations at each value of n to be tested is still required. 

For each simulation we make the statistic 
R = 



where Cf{6a) is the value of the ACF of a simulation at bin A and Cd{6a) is the value of 
the ACF of the data at correlation bin A. If R = 0, the ACF of the simulation is the same 
as that of the data. 

We calculate the likelihood, P(D\n, Q), or the relative probability of obtaining the 
FIRS data ACF given a model parametrized by n and Q, by choosing a limiting value, 
Rum, for R and finding the number of simulations, N, for which R < Ru m . We normalize 
the distribution so that the maximum likelihood is unity. The a posteriori probability that 
the Universe is parameterized by a certain n and Q can then be found by invoking Bayes's 
theorem with a suitable a priori probability distribution for (n, Q) (see, for example, Martin 
(1971)). 

In essence, we are treating the correlation functions as vectors in an n c dimensional 
ACF space. The various i?'s represent the weighted magnitude of the difference between 
the data ACF vector and a simulation ACF vector. If the T^'s were all the same, Ru m 
would define the radius of an n c dimensional sphere centered upon the data ACF. When 
the Ta are different, the sphere is deformed into an ellipsoid. By counting the number of 
simulations with R less than some chosen Rum, we are calculating the density of simulated 
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ACF's in the neighborhood of the measured ACF. This density is proportional to the 
likelihood. The method accounts for non-uniform sky coverage, instrumental noise, cosmic 
variance and the effects of removing an offset and dipole from an the data. 



Figure 2 shows the results of the analysis for — 1 < n < 3 and /iK < Q < 50 /iK. 
In this case, Ru m = 430 ( /iK) 2 , resulting in N max = 171 simulations at the most likely 
combination of (n, Q) with R < Ru m . That is, for (n, Q) = (1.0, 19 /iK), 171 out of a total 
of 20000 simulations had R less than 430 ( /iK) 2 . All other combinations of (n, Q) had fewer 
simulations fall within the test volume. Figure 2 shows likelihood contours for 0.05, 0.25, 
0.5 and 0.75 with solid lines along with a large x at the maximum. The inner broken line 
marks the P(D\n, Q) = e _1//2 likelihood contour, while the outer broken line corresponds to 
P(D\n, Q) = 0.34. They are shown for comparison to other analyses. If one were to extract 
only the n — 1 portion of Figure 2, in the limit that the extracted likelihood is Gaussian, 
the P(D\n, Q) = e -1 / 2 points would equal the la credible limits assuming a uniform prior. 
The P(D\n, Q) = 0.34 contour represents the 68% credible limits if we use a prior of one 
for the models considered here and zero for all others. In other words, the the sum of the 
likelihoods within the 0.34 contour is 68% of the sum of all the likelihoods in Figure 2. 

For each model, the number of simulations within our limiting radius is governed by 
the binomial distribution. The probability of a particular simulation falling within the 
test volume is much smaller than that of it falling outside the volume. We can, therefore, 
approximate the variance in the estimate of the various iV's by N itself. Hence, the 
uncertainty in the probability for those models with P(Q,n\D) 1 is approximately 



y2/N max , and the error in the estimates of the likelihoods near the maximum likelihood in 
Figure 2 is approximately 11% (with iV max = 171). This is borne out by Figure 3, which 
shows that the (n, Q) with the maximum likelihood depends slightly on the values of Ru m , 
but not by amounts in excess of what is expected statistically. 

Bond (1994) has noted that the FIRS data contain "white noise." In principle, this 
noise affects only the Q\ — bin of the correlation function, which is not included in this 
analysis. To check the possibility that the noise has 'leaked' out of the zero ACF bin, 
the analysis was repeated after excising the bin at 8a = l-°4 (effectively eliminating all 
correlations on angular scales less than 2.°8). The results are consistent with those quoted 
above. If the data contained only white noise, Figure 2 would have a maximum at n = 



4. 



Results and Discussion 
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3. Clearly, the white noise does not significantly contaminate these results. Again, this is 
because this method is primarily sensitive to spatial correlations in the data. 

This method has been checked in a number of ways. Figure 3 shows that there are 
small changes in the likelihood contours as Ru m is increased. This indicates that a smaller 
value of Rum is desirable. With this number of simulations, however, reducing Ru m degrades 
the plot to the point where it is no longer useful. Note though that the likelihood is stable 
for n ~ 1 and that for other n's the likelihood values are conservative (that is, the likelihood 
is overestimated). 

We have also checked that changing the weighting in the definition of R does not affect 
the results. Redefining R such that 



where Wa is now defined as 1, s/Ta, T\ or T\ yields consistent results, though, the bounds 

are not as strict. 

We put limits on the bias of this method by applying the procedure to 500 additional 
simulations with (n,Q) = (1.0,20 /iK). Though not a well defined quantity because some 
of the simulations would naturally prefer an n < —0.9 or an n > 2.9, the mean result, 
(n,Q) = (0.8, 16.5 fiK), is consistent with the input values to within the uncertainties in 
estimating the mean. 

The approximation of the FIRS beam by a 3.°8 Gaussian is adequate. Results 
from simulations made with a Legendre transform of the measured profile do not differ 
substantially. If a Gaussian of width 4.°2 is assumed to account for smearing, then the 
value for Q at n = 1 increases by 1 /iK. 

Figure 2 implies that for the FIRS data n « 1 is favored. We point out, however, that 
only a fraction of the total ACF space has been tested. Many other possibilities exist. 
Wright et al. (1992) found the ACF of the COBE/DMR first-year data to be well described 
by a Gaussian with a coherence angle, 6 C , of 13. °5; that is, the correlation function of the 
DMR data was fit well by a function of form C(9) = C(0)e~°' 5 ^/ e ' c ) 2 with adjustments for 
offset and dipole removal, where C(0) is the variance of the intrinsic sky fluctuation. This 
statistical method is well suited to testing this possibility. Retaining the normalization 
used in Figure 2, we find the likelihood of such a model with C(0) « (44 fiK) 2 to be 1.4. 
That is, 40% more likely than the most likely model based on a power law spectrum. Thus, 
while n 1 may be preferred compared to other models with a power spectrum of density 
perturbations that follows a power law, we cannot prove it is the best model. 
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We emphasize that no one experiment can prove that the correlated signal measured is 
due solely to the CMB. However, in the case of the FIRS data we have ruled out correlation 
stemming from astrophysical foregrounds, the instrument, and analysis artifacts (Page et al. 
1994b). The similarity of this result to the COBE/DMR auto-correlation function and 
the significant cross-correlation (Ganga et al. 1993) with DMR strongly suggest that both 
experiments are detecting correlations in the CMB. Work is underway to generalize our 
method to find the most likely (n,Q) for the FIRS/DMR cross-correlation. 
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6. Figure Captions 



Fig. I. — The FIRS auto-correlation function. The error bars represent the variations due 
to experimental noise only. For this plot, the 180° span of 9 a is divided into 40 bins each 
of width 4.°5. Due to the limited sky coverage, only 36 bins contain data. The shaded 
area shows the la spread in the 20000 simulated auto-correlation functions for n = 1 with 
an amplitude of Q = 19 fiK. The broken curve (corresponding to the right axis) shows the 
number of pairs of pixels contributing to the correlation function at each correlation function 
bin. 

Fig. 2. — Likelihood contours for Ru m = 430 ( fiK) 2 , corresponding to a maximum of 171 
simulations included within the test volume. The solid contours correspond to likelihoods of 
0.05, 0.25, 0.5, and 0.75. The x at (n, Q) (1.0, 19 fiK) is the maximum and the broken 
lines correspond to likelihoods of 0.34 and 0.68. 

Fig. 3. — Likelihood contours for various values of Ru m - The values of Ru m are 410, 420, 430 
and 440 ( fiK) 2 , moving clockwise from the top left. The maximum number of simulations 
within the test volumes for each case is indicated on the plots. 
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